Explosive percolation: a numerical analysis 
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Percolation is one of the most studied processes in statistical physics. A recent paper by Achliop- 
tas et al. [Science 323, 1453 (2009)] has shown that the percolation transition, which is usually 
continuous, becomes discontinuous ( "explosive" ) if links are added to the system according to special 
cooperative rules (Achlioptas processes). In this paper we present a detailed numerical analysis of 
Achlioptas processes with product rule on various systems, including lattices, random networks a la 
Erdos-Renyi and scale-free networks. In all cases we recover the explosive transition by Achlioptas 
et al.. However, the explosive percolation transition is kind of hybrid as, despite the discontinuity 
of the order parameter at the threshold, one observes traces of analytical behavior, like power law 
distributions of cluster sizes. In particular, for scale-free networks with degree exponent A < 3, all 
relevant percolation variables display power law scaling, just as in continuous second-order phase 
transitions. 

PACS numbers: 89. 75. He, 05.45.Df 
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I. INTRODUCTION 

Percolation phenomena [l[ represent probably the sim- 
plest examples of phase transitions that one could pos- 
sibly imagine. On infinite lattices, the process consists 
in occupying sites or bonds/links with some probability 
p. Nearest-neighboring occupied sites or links form clus- 
ters. When p exceeds a given system-dependent thresh- 
old value p Cl a macroscopic cluster, i.e. a cluster oc- 
cupying a finite fraction of all available sites or links, 
is formed {percolation cluster). This transition is con- 
tinuous, or second-order, as the order parameter varies 
smoothly from zero to values greater than zero. The same 
type of connectedness transition not only occurs on reg- 
ular graphs like lattices, but on any type of graphs. On 
random networks a la Erdos-Renyi (ER) 0] , for instance, 
one starts from a set of N nodes and adds links such that 
the probability p that two nodes are joined by a link is 
the same for all pairs of nodes. When p exceeds the value 
p c ~ a percolation cluster, or giant component, 

emerges and the transition is again continuous. Another 
well studied example is that of random networks with 
power law degree distributions of degree (number of node 
neighbors), usually called scale-free (SF) networks 
Here the process is better defined by removing, rather 
than adding, links. Links are removed until the graph 
is fragmented into microscopic clusters, i.e. there is no 
giant component. Remarkably, it has been shown that, 
if the exponent A of the degree distribution is smaller 
than 3, the giant component disappears only if one re- 
moves nearly all links of the graph, so that the fraction of 
remaining links with respect to the initial number goes 
to zero in the limit of infinite system size (6l-[l0|. This 
can be equivalently stated by saying that the percolation 
threshold is zero. Nevertheless, whether the threshold is 
zero or non-zero, the percolation transition is still contin- 
uous. In fact, the continuous character of the transition 
is a feature of all known percolation processes. However, 
this is true for random percolation, where links are ran- 



domly placed on the system, like in the examples above. 
Recently, Achlioptas et al. have shown that, if links are 
placed according to special cooperative rules, the perco- 
lation transition may become discontinuous [ll|. Such 
rules are non-local in character, as they require informa- 
tion between different parts of the system. Achlioptas et 
at introduced their rules in the growth of random net- 
works, and found an abrupt jump in the size of the giant 
component at the percolation threshold, hence the name 
explosive percolation. This peculiar type of transition is 
due to the fact that links are placed such to consider- 
ably slow down the formation of large clusters, so that 
clusters are mostly of about the same size [H, EH • In 
this way one reaches a point in which the insertion of a 
vanishingly small fraction of links leads to the merger of 
most of such small clusters, generating a big macroscopic 
cluster. The same effect has been observed by Ziff on 2- 
dimensional lattices [l|J|. For SF networks, the problem 
has been studied by Cho et al. [l5[ and by the authors 
of this paper [l6j], with different conclusions. In Ref. [l5[ 
the authors conclude that the percolation transition is 
discontinuous for any value of the degree distribution ex- 
ponent A greater than a critical value A c ~ 2.3; we found 
that, for A c < A < 3, the behaviour at the percolation 
threshold is consistent with that of a continuous transi- 
tion, while for A > 3 the expected behavior of a discon- 
tinuous transition is recovered. 

In this paper, we carry out an extensive numerical 
analysis of the phenomenon of explosive percolation. We 
will describe the case of SF networks [l6| , which we have 
studied in our previous paper, but we will also present 
results on lattices and random networks. The results of 
random percolation in all graph topologies will be pre- 
sented too, for comparison. Like in the paper by Achliop- 
tas et al., all graphs discussed in this paper will be built 
through dynamic growth processes. 

The paper is organized as follows. In Section|TTl we de- 
scribe the growth models considered in this paper. Sec- 
tion IHT] contains the results obtained from our numerical 
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Figure 1: (Color online) Scheme of an Achlioptas process with 
product rule. One of the two links represented by the dashed 
lines, which are selected at random among all possible pairs of 
non-adjacent nodes, has to be eventually added to the system. 
According to the product rule, the winner is the link that joins 
the pair of clusters with the smaller product size. In this case 
the winning link is that between clusters ci and C2, whose 
product size is 7 ■ 2 = 14, which is smaller than the product 
of the sizes of clusters C3 and C4 (4 • 4 = 16). 



simulations. In Section H*Vl we discuss the results of the 
analysis. A summary is presented in Section IVl 



II. GROWTH MODELS 

Our simulations of random percolation will be per- 
formed according to the Random Growth (RG) model, 
i.e., by iteratively adding one link to a system with N 
nodes, where the link is rando mly selected among all 
possible links. This procedure [lT], Gil is equivalent to 
classical bond percolation. 

Actually, also in Achlioptas processes links are added 
one by one. The difference is that the link to be added 
is chosen among two or more randomly selected links, 
according to a deterministic rule. In this paper we focus 
on the Product Rule (PR), which prescribes that the link 
to be picked is the one minimizing the product of the sizes 
of the two clusters joined by the link. The process is 
schematically illustrated in Fig. [TJ In the paper we shall 
often call this specific Achlioptas process as PR model, or 
simply PR. Other options are available too. For instance 
one could go for the sum of the cluster sizes, instead of the 
product. If, among two or more randomly selected links, 
the choice is random, one recovers random percolation. 
The systematic minimization criterion slows down the 
process of cluster growth, decelerating the percolation 
transition, which then may become "explosive" . 

For both RG and PR models, the growth proceeds until 



one reaches the desired density of links p. We defined p 
as the number of links of the graph divided by the total 
number of links present in the graph when it has been 
"completed", i.e., when the last link has been added. All 
graphs considered in this paper are "sparse", i.e., the 
ratio (k) between (twice) the number of links and the 
number of nodes N, which is the average degree of the 
nodes, does not depend on N. Therefore, since at time 
t of the growth process there are exactly t links in the 
system, their density p, according to our definition, is 
t/(N(k)). 



III. NUMERICAL ANALYSIS 

Our numerical analysis aims to understand and char- 
acterize the nature of the percolation transition induced 
by RG and PR. In order to do that, we make use of fi- 
nite size scaling [l9| |. a well-known technique adopted in 
numerical studies of phase transitions. For continuous 
phase transitions, every variable X near the threshold p c 
is scale- independent, due to the infinite correlation length 
of the system at p c , so it has a power law form, 



X ~ \p - p c 



(1) 



where w is a critical exponent. On a finite system of size 
N, the variable X has the following scaling form near the 
threshold 
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X = N 



(p- Pc ) N 1 '"' 



(2) 



where v is another critical exponent and F a universal 
function. For p — p c , the variable displays the simple 
scaling X ~ AT-"/", which can be used to deduce the 
exponents' ratio oj/v, from the examination of several 
systems with different sizes. Also, if p Cl v and u> are 
known, by plotting the expression X~N^I V as a function of 
(p — p c ) N 1 /" one yields the universal function F, which 
does not depend on N, so curves corresponding to differ- 
ent system sizes collapse. 

In this work we examined the two main variables mea- 
sured in percolation, i.e. the percolation strength P and 
the average cluster size S. The percolation strength P 
is the order parameter of the transition, and measures 
the relative size of the percolating cluster (s) with respect 
to the total system size N. On generic graphs there is 
no operative criterion to define a percolating cluster (as 
opposed to lattices), so one usually takes P as the rel- 
ative size of the largest connected cluster. The critical 
exponent of the percolation strength is indicated with (3 
and the scaling ansatz of P is 



p = N -P/v F (i) Up- Pc ) N 1 '" 
The average cluster size S is defined as 



S = 



(3) 



(4) 
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where n s stands for the number of clusters of size s per 
node. The sums run over all possible values of s except 
for the one of the largest cluster. The critical exponent 
of the average cluster size is indicated with 7 and the 
scaling ansatz of S is 



S = N l/v F (2) \{p-p c ) N 1 /" 



(5) 



The universal functions and of Eqs. ((3]) and ([5|) 
are different from each other, although they are related. 
We remark that in Ref. [TJ[ we have used the suscepti- 
bility x °f the order parameter, which measures the size 
of its fluctuations, rather than the average cluster size 
S. Therefore, the values of 7 that we present here are 
different from those of Ref. [16| ■ 

In random percolation, the probability distribution 
P{s) of sizes of the "finite" clusters, i.e. of all clusters ex- 
cept the largest, decreases as the power law P(s) ~ s~ T 
with the system size s at the percolation threshold. In 
our simulations we have also measured the critical expo- 
nent t (usually called Fisher exponent). We remark that, 
for a given system, P(s) is proportional to n s . Their re- 
lation is P(s) — Nn s /n Cl where n c is the total number 
of "finite" clusters. This is why in the paper we shall 
use the symbol n s to indicate P(s) as well. In the plots, 
however, n s is normalized as P{s), for consistency. 

In lattice percolation, as well as in spin models, the 
exponents /3l, and (where L stays for lattice) are 
linked by the so-called hyperscaling relation 
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where d is the dimension of the lattice. In the general 
case of graphs, we do not have a space dimension, so the 
scaling is done in terms of the "volume" N, as we have 
done in Eqs. (J2J), © and ([5]). In lattices N — L d and 
the hyperscaling relation for the exponents expressing the 
scaling of the variables in terms of the volume N reads 



7 , 2/3 



1. 



(7) 



Eq. ([7} is actually very general, and holds for random 
percolation on any system below the upper critical di- 
mension [20| . 

The identification of the percolation threshold p c is 
performed in two independent ways. One way consists in 
using the scaling of the pseudo-critical points p c (N) 



P, 



= Pc {N)+bN- 1 '\ 



(8) 



where b is a constant which has to be determined from 
the fit together with the other parameters v and p c . The 
pseudocritical points can be defined in several ways. We 
took the positions of the peaks of S at different system 
sizes N. 

The second method is based on Eq. By plotting 
the percolation strength P as a function of the system 



size N for a given value of p, the correct value of the 
percolation threshold can be determined by finding the 
value of p which yields the best power law fit. 

In [lH a new method for the determination of the na- 
ture of the transition has been proposed. The method 
consists in studying the behaviour of the width of the 
transition window as a function of the system size. As a 
measure of the width of the transition window we con- 
sidered the quantity Ap = pi—p\ 1 where pi is the lowest 
value of p for which P > 0.5 and p\ the lowest value 
of p for which P > 1/^/N. As we will see, the width 
of the transition window generally scales as a power law 
with the system size and its dependence from N can be 
therefore written as Ap ~ N~ a . Achlioptas et al. argued 
that, for continuous transitions, Ap should be indepen- 
dent of the system size (a = 0), whereas, if there is an 
explosive first-order transition, Ap should decrease with 
N (a > 0). Actually, in 2-dimensional lattices Ziff has 
found that a > even in the case of random percola- 
tion [TJJ • This is however due to the fact that in the par- 
ticular case of the lattice P2 is essentially coincident with 
the actual critical threshold of the system; therefore on 
the lattice one should take a value p% appreciably larger 
than p c (like the point at which P > 0.7, for instance). 
We stress that the choice of p\ and P2 is completely ar- 
bitrary, so the robustness of the exponent a needs to be 
tested. Therefore we also used another definition of Ap, 
namely Ap = p2— Px, where p2 is the lowest value of p for 
which P > 0.2. Also in this case, we can generally write 
Ap ~ N~ a . The robustness of the scaling of Ap would 
be indicated by the equality of the exponents a and a. 



A. Lattices 

We consider first the case of 2e?-lattices (square lat- 
tices) with periodic boundary conditions. The results ob- 
tained from our simulations with RG (see Figure [2]) con- 
firm the well-known fact that the transition is continuous. 
We also recover the correct critical exponents; the scaling 
is done in terms of the linear dimension L of the lattice, 
as it is customary in this case. PR on 2<i-lattices has been 
only recently studied by Ziff [JJJ , who has shown that the 
transition is explosive, like that observed by Achlioptas. 
In Figs. [2b, Hbl, and[5f we report the results obtained 
by applying PR on 2d-lattices. We find a trivial scaling 
for the order parameter P, with exponents' ratio f3 /v ba- 
sically equal to zero [0.07(3)] (Fig.[2t). This is consistent 
with what one expects to find for a discontinuous transi- 
tion. On the other hand, we find a clean non-trivial power 
law scaling at p c for the average cluster size S, with ex- 
ponent j/u = 1.7(1) (Fig. [5J1). This had been observed 
by Cho et al. in SF networks [l5j]. An explanation of 
this is provided by Fig. [5fe, which shows the distribution 
of sizes n s for all clusters except the largest one. The 
distribution is a clear power law [exponent 1.9(1)], which 
is unexpected for a classic discontinuous transition, as it 
usually occurs for continuous transitions. Therefore, all 
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variables derived from n s , like the average cluster size S, 
display power law scaling. This striking feature, as we 
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Figure 2: (Color online) Analysis of 2d-lattices. (a) RG 
model: the percolation strength P is plotted as a function 
of the lattice side L for three different values of the occupa- 
tion probability: p = 0.499 (violet diamonds), p — 0.5 (orange 
circles) and p = 0.501 (grey squares). The dashed line stands 
for the best fit obtained at the critical point p — p c = 0.5, 
which allows to determine f3/v = 0.11(1). (b) RG model: the 
average cluster size S is plotted as a function of the lattice 
side L for the same values of p as those used in (a). The 
dashed line has slope 7/1/ = 1.76(1). (c) PR model: the 
percolation strength P is plotted as a function of the lattice 
side L for three different values of the occupation probabil- 
ity: p — 0.5256 (violet diamonds), p = 0.5266 (orange circles) 
and p — 0.5276 (grey squares). The dashed line stands for 
the best fit obtained at the critical point p = p c — 0.5266(2), 
which allows to determine f3/v — 0.07(3). (d) PR model: S is 
plotted as a function of the lattice side L for the same values 
of p used in (c). The dashed line has slope 7/1/ = 1.7(1). (e) 
Comparison between the cluster size distributions measured 
at the critical threshold for both growth models. For RG (or- 
ange circles) r = 2.05(1) (black dashed line), while for PR 
(grey squares) r = 1-9(1) (red dotted line). Simulations have 
been performed on systems with L = 4096. (f) Ap as a func- 
tion of the system size N: a = 0.15(1) (dashed black line) for 
RG (orange circles) and a — 0.24(1) (dotted red line) for PR 
(grey squares). The first value is questionable, as the scaling 
should yield a plateau (a ~ 0), like we have indirectly veri- 
fied (see text). To see the correct scaling one should simulate 
much larger systems. 

will see below, is common to all "explosive" transitions 
we have investigated here. 

Finally, in Fig. [5f we show the results of the Achlioptas 



test for both RG and PR. For RG, we find a = 0.15(1). 
As we remarked above, a is non-zero despite the contin- 
uous percolation transition, which seems to go against 
the argument by Achlioptas et at. However, this hap- 
pens because P2 is very close to the critical point p c . The 
correct behavior can be seen if one considers a window 
clearly including p Cl which could be done by taking a 
larger value for the upper limit of the window, like, e.g. 
the smallest value p% at which the relative size of the giant 
component exceeds 0.7. Actually the scaling of pz — p\ 
(not shown) still shows sublinear behavior, but we believe 
that this is due to the fact that p\ grows too rapidly for 
the systems we were able to simulate. In fact, p% — pi 
is approximately constant for the lattice sizes we have 
taken, so p% — p\ > P3 — pi cannot go to zero in the limit 
of infinite lattice size. For PR, we obtain a = 0.24(1). 
This result is quite different from the value found by Ziff 
(0.34). However, in his simulations, Ziff has considered 
only links between clusters, whereas we have considered 
all possible links, including those within clusters. Sim- 
ulations of the process a la Ziff has confirmed that this 
is indeed the reason of the discrepancy with our result. 
We have performed the same analysis for the window Ap 
defined in Section HTT1 the exponents a for RG and PR 
are consistent with the corresponding values of a (see 
Table QJ. 

In 3d-lattices, the general picture is consistent with 
that in two dimensions (Fig. [3]). Classic percola- 
tion results, threshold and exponents, are recovered 
(Fig. [Sh., [3b, [5b). The scaling at p c of the order parame- 
ter P for the PR process is again trivial, with exponent 
P/v = 0.02(2), essentially zero (Fig. [3b). The S scales 
with an exponent 7/1/ = 2.1(1) (Fig. [5bl), again due to 
the power law shape of the distribution of cluster sizes 
(Fig. [3b^. We remark that the exponent r = 1.99(4) is 
compatible with that we found in two dimensions [1.9(1)]. 
The test of Achlioptas et al. (Fig. [5f) yields again a 
non-zero value of the exponent a for RG (a = 0.10(1)) 
(probably because our lattices are not yet large enough 
to see the actual behavior, as in 2c?), and a larger value 
for PR (a — 0.30(1)). Like in two dimensions, also in 3c? 
the exponents a for RG and PR are consistent with the 
corresponding values of a (see Table HJ. 



B. Erdos-Renyi random networks 

Percolation studies on random networks a la Erdos- 
Renyi (ER) have a long tradition, as we wrote in the 
Introduction. Fig. 3] summarizes the results of our anal- 
ysis. The well-known results of random percolation, 
threshold and exponents, are recovered, as illustrated in 
Figs.[3bo,[3}3 and0b- In particular, we notice that the hy- 
perscaling relation of Eq.[7|is satisfied for the exponents' 
ratios j3/f and j/v. For PR, instead, we see again a flat 
profile of the order parameter P with (fi/v — 0.02(1)), 
which hints to a discontinuous transition, together with a 
power law scaling of the average cluster size S, with expo- 
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Figure 3: (Color online) Analysis of 3d-lattices. (a) RG 
model: the percolation strength P is plotted as a function 
of the lattice side L for three different values of the occupa- 
tion probability: p = 0.2478 (violet diamonds), p = 0.2488 
(orange circles) and p = 0.2498 (grey squares). The dashed 
line stands for the best fit obtained at the critical point 
p = p c = 0.2488(3), which allows to determine fi/v = 0.48(1). 
(b) RG model: the average cluster size S is plotted as a func- 
tion of the lattice side L for the same values of p used in (a). 
The dashed line has slope y/v = 2.0(1). (c) PR model: the 
percolation strength P is plotted as a function of the lattice 
side L for three different values of the occupation probability: 
p = 0.3866 (violet diamonds), p — 0.3876 (orange circles) and 
p = 0.3886 (grey squares). The dashed line stands for the best 
fit obtained at the critical point p = p c — 0.3876(2), which al- 
lows to determine f3/v = 0.02(2). (d) PR model: the average 
cluster size S is plotted as a function of the lattice side L for 
the same values of p used in (c). The dashed line has slope 
y/v = 2.1(1). (e) Cluster size distributions n s for 3d-lattices 
at the critical point. For both percolation models n s ~ s~ T 
as s increases. For RG (orange circles) r = 2.20(1) (black 
dashed line), while for PR (grey squares) r = 1.99(4) (red 
dotted line). Simulations have been performed on systems 
with L = 256. (f) We plot the quantity Ap as a function 
of the system size N. For both models Ap decreases as a 
power law, Ap ~ N~ a , as N increases. In particular we have: 
a = 0.10(1) (dashed black line) for RG (orange circles) and 
a — 0.30(1) (dotted red line) for PR (grey squares). The first 
value is questionable, as the scaling should yield a plateau 
(a ~ 0), like we have indirectly verified (see text). To see the 
correct scaling one should simulate much larger systems. 



Figure 4: (Color online) Analysis of ER random networks, 
(a) RG model: the percolation strength P is plotted as a 
function of the system size N for three different values of the 
occupation probability: p — 0.495 (violet diamonds), p — 0.5 
(orange circles) and p = 0.505 (grey squares). The dashed 
line stands for the best fit obtained at the critical point p = 
p c = 0.5, which allows to determine fi/v — 0.33(1). (b) RG 
model: the average cluster size S is plotted as a function 
of the network size N for the same values of p used in (a). 
The dashed line has slope y/v = 0.34(1). (c) PR model: the 
percolation strength P is plotted as a function of the network 
size iV for three different values of the occupation probability: 
p = 0.8872 (violet diamonds), p — 0.8882 (orange circles) and 
p — 0.8892 (grey squares). The dashed line stands for the 
best fit obtained at the critical point p = p c = 0.8882(2), 
which allows to determine /3/u = 0.02(1). (d) PR model: the 
average cluster size S is plotted as a function of the network 
size iV for the same values of p used in (c). The dashed line 
has slope y/v = 0.48(4). (e) Cluster size distributions n s for 
ER random networks at critical point. For both percolation 
models n 3 ~ s~ T as s increases. For RG (orange circles) r = 
2.51(2) (black dashed line), while for PR (grey squares) r = 
2.08(5) (red dotted line). Simulations have been performed 
on systems with N = 8192. (f) We plot the quantity Ap as a 
function of the system size N. For both models Ap decreases 
as a power law, Ap ~ N~ a , as L increases. In particular we 
have: a — 0.03(1) (dashed black line) for RG (orange circles) 
and a = 0.36(1) (dotted red line) for PR (grey squares). 



nent y/v = 0.48(4). The exponent t = 2.08(5) (Fig. Si) 
is still compatible with the values found for PR on both 
2d and 3d lattices (see Section [TiTaI and Table HJ. The 
Achlioptas test of Fig. 0f yields a — 0.03(1) for RG, com- 
patible with a window Ap that is independent of N, while 
for PR a — 0.36(1), in agreement with the calculations 
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Figure 5: (Color online) Achlioptas process with PR on ran- 
dom SF networks. The plot shows the percolation threshold 
Pc{N) as a function of the degree exponent A for various net- 
work sizes N. The black line represents the infinite size limit 
extrapolation of the critical threshold, performed by apply- 
ing Eq. J8). The percolation threshold becomes non-zero for 
A > Ac ~ 2.3. Reprinted figure with permission fromRef. 0. 
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Figure 6: (Color online) Percolation transition induced by 
an Achlioptas process with PR on SF networks. The degree 
exponent A = 2.5. (a) The percolation strength P is plotted 
as a function of the system size N for three different values 
of the occupation probability: p = 0.0529 (violet diamonds), 
p — 0.0629 (orange circles) and p = 0.0729 (grey squares). 
The dashed line stands for the best fit obtained at the critical 
point p = p c = 0.0629(1), which allows to determine Pjv = 
0.59(1). (b) The average cluster size S is plotted as a function 
of the network size N for the same values of p used in (a). 
The dashed line has slope 7/f = 0.24(1). 
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Figure 7: (Color online) Percolation transition induced by 
an Achlioptas process with PR on SF networks. The degree 
exponent A = 2.8. (a) The percolation strength P is plotted 
as a function of the system size N for three different values 
of the occupation probability: p = 0.1229 (violet diamonds), 
p = 0.1329 (orange circles) and p = 0.1349 (grey squares). 
The dashed line stands for the best fit obtained at the critical 
point p = p c = 0.1329(1), which allows to determine P/v — 
0.50(1). (b) The average cluster size S is plotted as a function 
of the network size N for the same values of p used in (a). 
The dashed line has slope 7/1/ = 0.42(1). 
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of Achlioptas et al. [ll| . Again, the same test performed 
with the window Ap yields essentially the same values 
of the exponent for both RG and PR (Table HJ, so the 
results of the test appear to be quite robust. 



Figure 8: (Color online) Rescaling of the percolation vari- 
ables P and S near the percolation transition induced by an 
Achlioptas process with PR on SF networks with degree ex- 
ponent A = 2.8. System size N goes from 256000 to 4096000 
via successive doublings. 



C. Scale- free networks 

SF networks have been objects of intense investigations 
over the last few years [^-Q . The main reason of their 
success is that they are a proxy of many natural, social 
and man-made systems, if the latter are represented as 
graphs. The ubiquity of networks with skewed distri- 



butions of degree is not accidental. Such broad distri- 
butions indicate that there is a whole hierarchy of node 
roles based on their degrees, going from a large majority 
of nodes with low degree to a small subset of nodes with 
high degree, or "hubs". The hubs have a fundamental 
role for the structure and dynamics of networks. Ran- 
dom SF networks with degree exponent A < 3 have so 
many hubs that a very small fraction of links (vanishing 
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in the limit of infinite system size) is enough to keep a 
macroscopic fraction of nodes of the graph in the same 
connected component, which can be equivalently stated 
by saying that the percolation threshold is zero 
In Ref. 16] we have already studied Achlioptas processes 
with PR on random SF networks. Here we present some 
more detailed calculations and add significantly new ma- 
terial. 

The networks are constructed as follows. The start- 
ing point is a set of N nodes and given degree sequence 
{fci, fo, . . . , fcjv}- The degrees of the sequence are taken 
from a power law distribution with exponent A. We set 
the average degree (k) equal to 5. If links are placed ran- 
domly, the procedure can be carried out with the config- 
uration model (2lj . i.e. by connecting randomly selected 
pairs of stubs adjacent to the nodes, until no more stubs 
are available. This is actually the procedure we have 
adopted for the RG model. For PR, instead, at each it- 
eration we pick two pairs of stubs and apply the PR to 
find which pair of stubs has to be eventually joined in a 
link (the PR applies as we have schematically illustrated 
in Fig. HI). 

In the degree exponent's range A < 3 we will present 
only results referring to PR, due to the absence of a per- 
colation threshold for RG. A remarkable result found in- 
dependently in Refs. [l5[ and [lj| is that the percolation 
transition of the Achlioptas process with PR has a non- 
vanishing threshold for A > A c ~ 2.3 (Fig. [5]). 

In Figs. [6] and [7] we show the scalings at p c of P and S 
for A = 2.5 and 2.8, respectively. At variance with what 
we have seen in Sections MI Al and MI B[ here the scaling 
of P at p c is non-trivial, as P decreases with A as a 
power law in both cases. This appears inconsistent with 
the typical scenario of a discontinuous transition, which 
generally yields the trivial scaling we have observed in 
Figs.[2p, [3fc and St. We shall come back to this issue in 
Section HVl A clean power law scaling at p c is also found 
for S (Figs. and[7}}), although we have seen that the 
same happens for explosive discontinuous transitions as 
well. 

In Fig. |5] we show the rescaling of the variables P and 
S. The data collapses observed in Figs. [Sfc and |8f show 
the profiles of the universal scaling functions F^- 1 ' and 
F (2) of Eqs. © and ©. The results of the Achlioptas 
tests for A = 2.5 and 2.8 are shown in Fig. [SJ In each 
case we present the scaling of both Ap and Ap, to check 
for the robustness of the results. We find that, while 
the scaling is clear for both variables, a =/= a. In fact, 
the exponents often indicate contradictory trends, with 
the transition window increasing and decreasing with N, 
which is clearly inconsistent. We cannot exclude that this 
is due to finite size effects (as we have seen on lattices) 
and that simulations on much larger systems would show 
consistent results instead. On the other hand, it might 
be that the results of Achlioptas test indeed depend on 
the specific definition of the transition window. 

For A > 3, however, the situation is different. Fig. [TU] 
reports the results of our finite size scaling analysis for 
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Figure 9: (Color online) Achlioptas test for SF networks, (a) 
A = 2.5: we plot the quantity Ap as a function of the system 
size N. For both percolation models Ap decreases as a power 
law, Ap ~ N~ a , as N increases. In particular we have: a = 
—0.04(1) (dashed black line) for RG (orange circles) and a = 
—0.26(3) (dotted red line) for PR (grey squares). We also 
consider the transition window Ap (defined in Section |III|) . 
from which we obtain: a = —0.06(1) (lower black dashed line) 
for RG (orange diamonds); a = 0.31(1) (lower red dotted line) 
for PR (grey triangles), (b) A = 2.8: same plot as the one 
of (a). The measured exponents are: a = —0.04(1) (dashed 
black line) for RG (orange circles); a = 0.04(1) (dotted red 
line) for PR (grey squares); a — —0.07(1) (lower black dashed 
line) for RG (orange diamonds); a — 0.32(1) (lower red dotted 
line) for PR (grey triangles). 



percolation transitions induced by RG and PR on SF 
networks with exponent A = 3.5. In this case we also 
present the results of RG, because for A > 3 there is a 
non-zero threshold. In Ref. [2(| it has been proved that, 
for random percolation on SF networks, f3/v = 1/(A — 1) 
and j/p = (A - 3)/(A - 1) for 3 < A < 4. For A > 4, 
the process reaches the mean field limit and the expo- 
nents are frozen: /3/v — 7/1/ = 1/3. Interestingly, these 
are just the values of the exponents for the percolation 
transition of ER random networks. SF networks tend 
to ER random networks in the limit A — > 00. Our es- 
timates of the exponents' ratios @/v and j/u (Figs. [TUb , 
and llOb ) are consistent with the predicted values for ran- 
dom percolation presented above. For PR, instead, we 
recover the same scenario as on lattices and ER random 
networks. The scaling of P at p c is trivial (Fig. [TUbl. 
with f3/v = —0.06(3), which is essentially zero, while the 
power law scaling of S at p c is non-trivial (Fig. [TOHh 
with 7/1/ = 0.40(9). The Fisher exponent r = 2.2(1) 
(Fig. D]3i) . The Achlioptas test (Fig. [TUf) yields essen- 
tially the same set of values we had found for ER random 
networks (see Table JJ. Moreover, the values are stable 
no matter whether one uses Ap or Ap. 



IV. DISCUSSION 

In this section we want to discuss the results we have 
obtained, which are summarized in Table [U We have 
seen that our finite size scaling analysis leads to two dif- 
ferent scenarios. The first scenario is consistent with the 
"explosive" transition observed by Achlioptas et al., and 
occurs on ER random networks, lattices and SF networks 
with degree exponent A > 3. In all these cases we have 



8 



System 


Growth Model 


Pc 


P/v 


■y/u 


T 


Q 


a 


2d-lattice 


RG 


0.5 


0.11(1) 


1.76(1) 


2.05(1) 


0.15(1) 


0.16(1) 


PR 


0.5266(2) 


0.07(3) 


1.7(1) 


1.9(1) 


0.24(1) 


0.23(1) 


3d-lattice 


RG 


0.2488(3) 


0.48(1) 


2.0(1) 


2.20(1) 


0.10(1) 


0.10(1) 


PR 


0.3876(2) 


0.02(2) 


2.1(1) 


1.99(4) 


0.30(1) 


0.31(1) 


ER network 


RG 


0.5 


0.33(1) 


0.34(1) 


2.51(2) 


0.03(1) 


0.04(1) 


PR 


0.8882(2) 


0.02(1) 


0.48(4) 


2.08(5) 


0.36(1) 


0.36(1) 


SF network A - 2.5 


RG 





— 


— 


— 


-0.04(1) 


-0.06(1) 


PR 


0.0629(1) 


0.59(1) 


0.24(1) 


2.15(2) 


-0.26(3) 


0.31(1) 


SF network A - 2.8 


RG 











-0.04(1) 


-0.07(1) 


PR 


0.1329(1) 


0.50(1) 


0.42(1) 


2.13(6) 


0.04(1) 


0.32(1) 


SF network A - 3.5 


RG 


0.078(1) 


0.38(1) 


0.15(2) 


2.94(1) 


-0.02(1) 


-0.02(1) 


PR 


0.2224(2) 


-0.06(3) 


0.40(9) 


2.2(1) 


0.34(1) 


0.35(1) 



Table I: The table summarizes the results obtained from our numerical analysis. Percolation threshold and critical exponents 
are reported for each system and growth model analyzed. 



derived the same picture from finite size scaling, in par- 
ticular the saturation of the order parameter P at p c with 
the size of the system N. On SF networks with A < 3 
the situation looks different, as there we have observed 
a clear power law scaling of P at p c , just as one would 
expect to find in continuous transitions. Moreover, the 
pseudo-critical points also show the clean power law scal- 
ing of Eq. [8] for A < 3 (that is how the critical thresh- 
olds of Fig. O were derived), which usually happens for 
continuous transitions. This appears to contradict the 
conclusion of Cho et al., who claim that the transition 
is always discontinuous on SF networks [psl ]. Cho et al. 
have adopted the model by Chung and Lu [22| to build 
their networks, which is different from the procedure we 
used, but we have verified that the results obtained in 
this way are consistent with ours. 

However, the seemingly continuous transition we ob- 
serve for SF networks with A < 3 has the surprising and 
somewhat disturbing feature that the hyperscaling rela- 
tion of Eq. [7] is violated, as one can easily verify through 
the values of (3/u and 7/1/ reported in Table|T]for A = 2.5 
and A = 2.8. Such violation could imply that the transi- 
tion is not continuous after all. In order to test this, we 
have computed the distribution of the values of the order 
parameter P at the pseudo-critical point p c (N), for PR 
on SF networks with A = 2.1, 2.5, 2.8, 3.5. The results are 
reported in Fig. 111! The two horizontal bands visible in 
the main frame of each of the four panels indicate that P 
oscillates between two values at the pseudo-critical point, 
which means that the transition is discontinuous. Inter- 
estingly, this is also found for A = 2.1 < A c . We could 
not carry out the finite size scaling analysis in this case, 
because the percolation threshold vanishes in the infinite 
size limit, but the result on finite systems, as shown in 
Fig. [Tib , is the same as those for A > A c . We conclude 
that the percolation transition for an Achlioptas process 
with PR is discontinuous on SF networks, for any value 
of the degree exponent A. Therefore, based on the results 



of this analysis, we have to partially modify the conclu- 
sion we had drawn in Ref. [16j . where we had stated that, 
for A < 3, the transition is continuous. There actually is 
a discontinuous jump of the order parameter at p c : nev- 
ertheless, all relevant percolation variables display power 
law scaling at the percolation threshold for A < 3, in par- 
ticular Eqs. fl3J), dU) and © hold, just like in standard 
continuous transitions. Therefore, we hesitate to state 
that the transition is first- or second-order, as it looks 
like an unusual mixture of both. Therefore, the regime 
of SF networks for A < 3 is very intriguing and deserves 
further investigations. 

Furthermore, the explosive transition observed in the 
other cases, including the original transition discovered 
by Achlioptas et al, is not a standard discontinuous tran- 
sition neither. The most striking feature here is that 
the size distribution n s of the "finite" clusters at p c is a 
power law, not exponential or Gaussian as one usually 
observes in first-order phase transitions. This fact has 
the consequence that all variables computed by means 
of n s also display non-trivial power law scaling at p c , 
as we have seen with the average cluster size S. Inter- 
estingly, the Fisher exponent r for every transition we 
have investigated is very close to 2, and consistent with 
this value within errors. In Fig. [T^] we plot all distri- 
butions n s we have computed. Indeed, we see that the 
curves are strongly overlapping, and that only the curves 
corresponding to the anomalous discontinuous transition 
found in SF networks with A < 3 perhaps deviate from 
the general pattern, though very little. 

Another striking feature of our findings is the existence 
of a non-zero percolation threshold for SF networks for 
A > A c ~ 2.3 (Fig. [5]), in contrast with the fact that the 
threshold for random percolation is zero until A = 3. In 
Ref. [ll| , Cho et al. suggested an interesting explanation 
for this result. They noticed that, since in Achlioptas 
processes the networks are not constructed through the 
random addition of links, the degree distribution of the 
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Figure 10: (Color online) Analysis of SF networks with degree 
exponent A = 3.5. (a) RG model: the percolation strength 
P is plotted as a function of the system size TV for three dif- 
ferent values of the occupation probability: p — 0.074 (violet 
diamonds), p — 0.078 (orange circles) and p = 0.082 (grey 
squares). The dashed line stands for the best fit obtained at 
the critical point p — p c = 0.078(1), which allows to deter- 
mine )3/v — 0.38(1). (b) RG model: the average cluster size S 
is plotted as a function of the lattice side L for the same val- 
ues of p used in (a). The dashed line has slope 7/1/ = 0.15(2). 
(c) PR model: the percolation strength P is plotted as a 
function of the system size TV for three different values of 
the occupation probability: p — 0.2214 (violet diamonds), 
p — 0.2224 (orange circles) and p = 0.2234 (grey squares). 
The dashed line stands for the best fit obtained at the critical 
point p = p c = 0.2224(2), which yields f3/v = -0.06(3). (d) 
PR model: the average cluster size S is plotted as a func- 
tion of the network size TV for the same values of p used in 
(c). The dashed line has slope 7/1/ = 0.40(9). (e) For both 
growth models n s ~ s~ T as s increases. For RG (orange cir- 
cles) p c = 0.078(1) and r = 2.94(1) (black dashed line), while 
for PR (grey squares) p c = 0.2224(2) and r = 2.2(1) (red dot- 
ted line). Simulations have been performed on systems with 
TV = 8192. (f) We plot the quantity Ap as a function of the 
system size TV. For both growth models Ap decreases as a 
power law, Ap ~ N~ a , as TV increases. In particular we have: 
a = —0.02(1) (upper dashed black line) for RG (orange cir- 
cles) and a = 0.34(1) (dotted red line) for PR (grey squares). 
We also consider the transition window Ap = p2—pi, where P2 
is the minimal value of the occupation probability at which 
P = 0.2. In this case we find again a good power law fit. 
For RG (orange diamonds) the decay exponent is unchanged, 
since a = —0.02(1) (lower black dashed line). Similarly for 
PR (grey triangles) a = 0.35(1). 



system during the growth deviates from that imposed by 
construction, which will be eventually reached at the end 
of the process. In Fig. Q2] we plot the degree exponent 
X e ff measured at the critical point as a function of the 
imposed exponent A. We see that the two exponents are 
quite different, and that there is a simple linear relation 
between them. In particular, we notice that X e ff ~ 3 
when A = A c ~ 2.3. Therefore, at the percolation thresh- 
old, SF networks constructed with an Achlioptas process 
with PR for A > A c are actually SF networks with degree 
exponent bigger than 3. For SF networks with degree ex- 
ponent bigger than 3, random percolation has a non-zero 
threshold, and this could be the reason of the non-zero 
threshold we observe for A > A c . However, we stress that 
the "effective" SF networks produced by Achlioptas pro- 
cesses are not random, so there is no a priori guarantee 
that one finds the same results as for random SF networks 
as the argument above implies. Still, one could expect 
some qualitative agreement. This is confirmed by the fact 
that X e ff ~ 4 when A = 3. For degree exponents larger 
than 4, random SF networks are hardly distinguishable 
from ER random networks. From the point of view of 
percolation the two classes of systems are in fact fully 
equivalent (same exponents). This could explain why, 
for A > 3, the picture we recover from finite size scaling 
looks the same as for ER random networks, and that the 
corresponding critical exponents (including the exponent 
a of the Achioptas test) are consistent with each other 
within errors (see Table [I]). 



V. SUMMARY 

In this paper we have performed a thorough numer- 
ical analysis of the percolation transitions induced by 
Achlioptas processes with product rule. The typical out- 
come, on lattices, ER random networks and SF networks 
is the "explosive" percolation transition originally ob- 
served by Achlioptas et al. [ll| . This transition is kind of 
hybrid, as it combines the discontinuity of the order pa- 
rameter at the critical point with analytical features like 
the power law decay of the size distribution of finite clus- 
ters, a feature typical of continuous transitions. Hybrid 
phase transitions are actually not new, they have been 
observed in a variety of domains, like spin glasses [23l. Izij . 
constraint satisfaction problems (K-SAT) |25l| and mod- 
els of jamming in granular materials [2614281. A remark- 
able feature of our findings is that the value of the ex- 
ponent r of the cluster size distribution appears to be 
compatible with 2 in all instances, despite the diversity 
of the systems we considered. 

For SF networks with degree exponent A < 3 the situ- 
ation is even more extreme: on the one hand, all perco- 
lation variables display power law scaling at the critical 
point, just like one expects for a continuous second-order 
phase transition; on the other hand, the order parame- 
ter still undergoes a discontinuous jump at the critical 
point. This is certainly something worth investigation in 
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Figure 11: (Color online) Achlioptas process with PR on SF networks. Distributions of the values of the order parameter P 
at the pseudocritical point p c (N) for different degree exponents A: 2.1 (a), 2.5 (b), 2.8 (c), 3.5 (d). The main frame of each 
plot shows the values of P for each of 1000 realizations. The insets display the distribution of the P-values (upper panel) 
and its cumulative (lower panel). The distributions are all bimodal, which indicates that the order parameter undergoes a 
discontinuous jump at the critical point. The network size is N = 8192000 in all cases. 



the future. As usually in numerical studies of phase tran- 
sitions, despite the large graph sizes we have investigated 
here, we cannot exclude that the regime we have tested 
is not yet "asymptotic" and that therefore is dominated 
by finite size effects, which give a distorted perception of 
what truly happens. We tend to discard this hypothe- 
sis, though, due to the remarkably clean scaling plots we 
have derived. 

Some theoretical arguments have be en prop osed to de- 
scribe explosive percolation transitions [12l . ll3l |. However, 
a real theory of such processes is still missing, and look- 
ing for a theory is certainly a challenging but promising 



future research direction. We hope that the results of our 
analysis will contribute to inspire new theoretical devel- 
opments in this topic. 
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